Quantum dynamics of one and two bosonic atoms in a combined tight-binding 

periodic and weak parabolic potential 
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Strongly interacting bosonic particles in a tight-binding periodic potential superimposed by a 
weak parabolic trap is a paradigm for many cold atom experiments. Here, after revisiting the sin- 
gle particle problem, we study interaction-bound dimers of bosonic atoms in the combined lattice 
and parabolic potential. We consider both repulsively- and attractively-bound dimers and find pro- 
nounced differences in their behaviour. We identify conditions under which attractive and repulsive 
dimers exhibit analogous dynamics. Our studies reveal that coherent transport and periodic os- 
cillations of appropriately prepared one- and two-atom wavepackets can be achieved, which may 
facilitate information transfer in optical lattice based quantum computation schemes. 

PACS numbers: 03.75.Lm, 37.10.Jk, 03.65.Ge 



I. INTRODUCTION 

Quantum transport in periodic structures is one of 
the central topics of condensed matter physics [1]. Re- 
cent interest towards spatially-periodic systems has been 
largely motivated by the remarkable progress in cooling 
and trapping bosonic and fermionic atoms in optical lat- 
tices 0- The relevant parameters of these systems can 
be controlled with very high precision and can be tuned 
to implement some of the fundamental models of con- 
densed matter physics. In a tight-binding regime, the 
Hubbard model accurately describes static and dynamic 
properties of these systems 0, Q . Importantly, in real ex- 
periments with cold atoms, the lattice is finite and often 
it is superimposed by a weak harmonic trap. This breaks 
the translational invariance of the lattice and thereby 
strongly modifies the properties of the system even in 
the limit of non- interacting particles [1, @, @] . In particu- 
lar, the low energy states of the single-particle spectrum 
behave like harmonic oscillator states while the higher en- 
ergy states are localized at the sides of the parabolic trap, 
which can lead to inhibition of quantum transport and 
dipole oscillations in a degenerate atomic gas [a 0, HI, @ • 

Atom-atom interactions profoundly enrich the Hub- 
bard model, as attested by, e.g., theoretical prediction 
followed by spectacular experimental demonstration of 
the transition from the superfluid to the Mott insulator 
phase in an ensemble of cold bosonic atoms in optical lat- 
tice [l(|. Strongly interacting bosons in periodic poten- 
tials can form tightly bound "dimers" [11| observed in a 
recent experiment |12j with repulsively interacting atoms 
in a lattice. Here we first discuss static and dynamic 
properties of a single atom in a combined periodic and 
weak parabolic potential and show that coherent trans- 
port and perfectly periodic oscillations of appropriately 
prepared wavepackets can be achieved. We then study 
the properties of strongly interacting atom pairs examin- 
ing both regimes of attractive as well as repulsive inter- 
actions. We show that, quite generally, the interaction- 
bound dimers behave as single particles with appro- 
priately rescaled parameters of the system. We iden- 



tify, however, important differences between attractively- 
bound and repulsively-bound dimers and find, rather sur- 
prisingly, that, in a weak trap, the repulsive dimer is 
bound stronger than the attractive one. 



II. THE MODEL 

We consider cold bosonic atoms in a combined tight- 
binding periodic and weak parabolic potential. In ID, the 
system is described by the Bose-Hubbard Hamiltonian 
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where fcj (bj) is the creation (annihilation) operator and 

hj = bjbj the number operator for site j, J is the 
tunnel coupling between adjacent sites, U is the on- 
site interaction, and £1 > quantifies the strength of 
the superimposed parabolic potential due to which site 
j = ±1,±2, . . . acquires energy offset fij 2 with respect 
to site j = corresponding to the minimum of the po- 
tential. A natural basis for Hamiltonian is that of 
the eigenstates \rij) = ~^( b j) n |0) of operator fij whose 
eigenvalues n = 0, 1, 2, . . . denote the number of particles 
at site j, and |0) = |{0j}) is the vacuum state. 

Single particle spectrum. We first discuss the single- 
particle case with the on-site interaction U playing no 
role. Recall that in the absence of parabolic potential, 
f2 = 0, the eigenstates of the Hubbard Hamiltonian ((T|) 
form a Bloch band of width 4J centered around zero. 
More quantitatively, given a finite flat lattice of N sites, 
the Bloch eigenenergies and corresponding eigenstates 
are given by 
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FIG. 1: Single particle energy eigenvalues Ek in a combined 
periodic and parabolic potential obtained by numerical diag- 
onalization of Hamiltonian |T} with J/Sl — 140. For compar- 
ison, the blue dots represent eigenvalues Ek within the Bloch 
(mini)band for a flat lattice (SI — 0) of length N — N — 31. 



with < k < N and Af a normalization constant. In 
the limit of N — > oo, Eq. (|2ap yields the well-known dis- 
persion relation [l[ E q = — 2Jcos(q) with < q < it 
and lattice constant d = 1. Remarkably, however, even 
a very weak parabolic potential !1« J drastically mod- 
ifies the spectrum of Hamiltonian JT]) 0, H[ , as shown in 
Fig. [1] Note that the spectrum, bound from below by 
—2 J, is composed of discrete energy levels Ek- Two dis- 
tinct groups of levels can be identified: (i) the low-energy 
levels Ek < 2 J forming a modified Bloch band, and (ii) 
the high-energy ones Ek > 2 J. 

(i) The parabolic potential effectively restricts the 
number of sites accessible to a particle with energy within 
the Bloch band — 2J < Ek < 2J. Roughly, only sites 
j = 0, ±1, for which Slj 2 < 2 J, can participate in 
the formation of the low-energy part of the spectrum || . 
More precisely, using second order perturbative correc- 
tions, we find that the modified Bloch band is restricted 
to sites j satisfying 

\j\<j^ ^1 + ^)^1.3^. (3) 

The low-energy part of the spectrum therefore contains 
N = 2 Umax J + 1 energy levels E , E X: . . .E N -x, which is 
illustrated in Fig. [H and we have verified this conclusion 
for a wide range of values of J/Vt. Note that the weak 
parabolic potential modifies the Bloch band in such a way 
that its lowest energy part is approximately linear in fc, 
similar to the spectrum of a harmonic oscillator. Indeed, 
using the properties of the Mathieu functions [HI , it can 
be shown [5| that in the limit of J /SI ^> 1, the low-energy 
eigenvalues Ek and eigenstates \xk) of Hamiltonian ([1]) 
are well approximated by 

E k w -2 J + 2-\/jQ (4a) 
\ Xk ) « AA^(2 fe fc!)- 1 / 2 e-^ 2 i/ fc (0)|l J ), (4b) 
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where AT is a normalization constant, Q = j ySY/J is 
the discretized coordinate, and Hk(C) is the fc-th Her- 
mite polynomial. Thus, the lowest energy eigenvalues Ek 
and eigenstates \xk) correspond to those of an energy- 
shifted harmonic oscillator with an effective frequency 
Huj = 2V JSl and an effective mass given by the usual 
expression /i = h 2 /(2JcP) [l[ which is valid near the bot- 
tom of the lowest Bloch band of a periodic potent ial. In 
particular, the ground state with Eq = —2J + V JSl is 
given by 

j 

The spectrum at the top of the modified Bloch band ap- 
proaches that of a uniform (SI — 0) finite (N = N) lattice 
given by Eq. (|2"a|) . 

(ii) On the other hand, the high-energy eigenvalues 
Ek > 2J with k > N are two-fold degenerate, as seen 
in Fig. [TJ The pairs of degenerate states with indices 
k = 2\j\ + 1 and k' — 2\j\ + 2 are localised around 
sites j = ±\j\ (\j\ > j max ) equidistant from the center of 
the parabolic potential, the corresponding energies being 
given by Ek.k< ~ Oj' 2 For such states, the localiza- 

tion occurs because, for large enough the transitions 
|lj) — » |lj±i) effected by the last term of Hamiltonian ([TJ) 
become non-resonant and the particle tunneling between 
neighboring lattice sites is suppressed. These high-energy 
states were shown to be responsible for damping of dipole 
oscillations and quantum transport of degenerate atomic 
gases in combined harmonic and optical lattice potentials 
[E B S B @- Such localised states can be selectively 
addressed by radio-frequency fields [l4[ and may be em- 
ployed for efficient initialization of a qubit register with 
fermionic atoms [15[ . 

Coherent dynamics of a single particle wavepacket. 
From the above analysis, it is clear that if we restrict 
ourselves to the harmonic oscillator-like states belong- 
ing to the lowest part of the energy spectrum, we can 
expect a quasi-periodic dynamics in the system. Non- 
dispersive transport of a single particle wavepacket from 
one side of the shallow parabolic potential to the other 
can then be achieved. In Fig.[2Ua) we show the dynamics 
of a single particle wavepacket represented by the 
ground state of the system |x ), Eq. ©, initially shifted 
by 7 sites from the trap center. Numerical solution of 
the Schrodinger equation using Hamiltonian ([1]) reveals 
almost perfect periodic oscillations of the discrete Gaus- 
sian wavepacket between the two sides of the parabolic 
potential with period r ~ 27r/w = (nh/J)y/J/Sl. 

From the set of harmonic oscillator-like states \xk) of 
Eq. (|4b[) . we can construct a well-localized wavepacket 
\ipV)) centered at a prescribed site j' (\j'\ < j max )- If we 
write the initial state as 

!#))> = $> fe |x fe >, (6) 

k 

the probability amplitude a.j for a particle to be at site j 
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FIG. 2: Time evolution of density pj = (hj) for a single parti- 
cle wavepacket \ip) in a combined periodic and parabolic po- 
tential with J/fi = 140. (a) Initial state |-0(O)) corresponds 
to the ground state |xo) (discrete Gaussian) shifted by 7 sites 
from the trap center, (b) Initial state \ip^ '(0)) is a localized 
around j' — 7 wavepacket constructed from the k = 0, 1, . . . 20 
eigenstates |xfc}- Insets in (a) and (b) show the density dis- 
tribution pj at t = and t ~ r/2. Time is measured in units 
of hj- 1 . 



is given by 

a, = (1^(0)) (x^A fe (2 fc fc!)- 1 /2 e -C J 2 / 2 i/ fe (C J ). (7) 
fe 

To obtain a localized around site j' state \ipv '), we max- 
imize \a,ji | 2 , which determines the set of coefficients {A^} 
in Eq. ©. In Fig.[3{b) we show the time evolution of such 
a localized state, which exhibits periodic collapses and 
partial revivals at sites —j' and j' with time steps r/2. 
The revivals are not complete since, as noticed above, the 
energy spectrum Ek for small k is only approximately lin- 
ear in k. Nevertheless, our results suggest that coherent 
non-dispersive transport of carefully engineered atomic 
wavepackets can be achieved in optical lattices in the 
presence of a shallow parabolic potential. 

III. TWO PARTICLE DYNAMICS 

We now consider two bosonic particles in the combined 
periodic and weak parabolic potential. Clearly, in the 
simplest case of feeble interaction \U\ -C J, we have two 
independent particles for which the results of the previ- 
ous section apply. But even for strong on-site interaction 
U , some aspects of the combined dynamics of two low- 




FIG. 3: Time evolution of density pj = (hj) for two particles 
in a combined periodic and parabolic potential with J/Q. = 
140 and U = — 10 J. (a) Initial state 1^(0)) corresponds to 
one particle in the ground state \xo) and the other particle in 
state |xo) shifted from the trap center by 7 sites, (b) Initially 
both particles in state |%o) are shifted from the trap center 
by 7 sites in opposite directions, (c) Initial state corresponds 
to both particles in state |xo) shifted from the trap center by 
7 sites in the same direction. Inset in (c) shows the projection 



energy particles can be inferred from the independent 
particle picture modified by short-range collisions. This 
applies when the initial state |^(0)) = \tp) <g> \ip') is com- 
posed of two non-overlapping single-particle wavepackets, 
| (ip | 2 -C 1, which upon collision with each other are 
reflected by the potential barrier \U\ > J. Examples of 
such a situation with large on-site attractive interaction 
energy U = —10J are shown in Figs.[3Ja) and[3Jb). Anal- 
ogous dynamics is observed for the repulsive interaction 
U = 10J. 

More intriguing is the case of initial state |^(0)) = 
\ip) <8> \tp) consisting of two overlapping single-particle 
wavepackets shown in Fig. G3c). This state has a signif- 
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icant population of the two-particle states \2j) given by 
J2j l(2j|^')| 2 — J2j l a j| 4 ' where <Zj are the single-particle 
probability amplitudes. Clearly, the population of two- 
particle states is largest in the central part of the initial 
density distribution. As seen in Fig. 02c), this part ex- 
hibits slow dynamics, characterized by the effective tun- 
nelling constant — —2J 2 /U (see below), and sepa- 
rates from the wings of the initial density profile. The 
wings, formed by the single-particle states oscil- 
late between the two sides of parabolic potential with 
the usual period r. 

Interaction-bound dimers. At this point, let us recall 
[Til OjJ that two bosonic particles occupying the same 
site j can form an effective "dimer" bound by the on- 
site interaction U. Thus, when \U\ 3> J, the first-order 
transitions \2j) — > |lj±i) effected by the last term of 
Hamiltonian (TT]) are non-resonant and the particles can 
not separate. However, the second-order in J transitions 
\2j) — > |2j±i) via virtual intermediate states |lj) |lj±i) 
are resonant. Consequently, the dimer can tunnel as a 
whole with the effective rate J (2) = -2J 2 /U « J [HI- 
This explains the dynamics seen in Fig. [3jc) where the 
initial density distribution splits into slow and fast prop- 
agating components, the former composed of the dimer 
states \2j) while the latter containing the monomer states 

1 1,')- 

If the initial state is prepared in such a way that 
only two-particle (dimer) states are populated, as im- 
plemented in, e.g., fl2| . for \U\ 3> J the system can, to a 
good approximation, be described by an effective dimer 
Hamiltonian derived in the second order in J/U fill ]. I n 
terms of the dimer creation = (b^) 2 [1 / \/2(fij + 1)] 
and annihilation Cj — [l/y / 2(hj + l)](bj) 2 operators, and 
number operator rhj — CjCj — fijj2, the effective Hamil- 
tonian for a single dimer reads [16[ 

H cS = [n^fm 3 + (U- J^)m 3 - JW(c]c 3+1 +c} +iCj )] , 
j 

(8) 

where — 2f2 is the strength of a parabolic poten- 
tial seen by the dimer, while (U — J^ 2 ') represents the 
"internal" energy of a dimer. 

Before proceeding, let us note that, differently from 
the flat lattice situation considered in [ll[, here the ef- 
fective Hamiltonian H c g is not applicable in the vicinity 
of sites \j\ ~ |[/|/(2J1) where near-resonant dissociation 
of a dimer can occur via transitions \2j) — > |lj) |lj±i). 
But since we are interested in the dynamics of low-energy 
dimers with \U\ ^> J ^> fi, such high-j states cannot be 
reached. 

Consider first the case of strong attractive interaction 
U < leading to a positive tunnelling constant > 0. 
Then the effective Hamiltonian ([8]) has the same form 
as the Hubbard Hamiltonian fTJ) for a single particle in 
a combined periodic and parabolic potential. We can 
therefore immediately write the lowest energy eigenvalues 
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FIG. 4: Time evolution of atom density pj = (hj) for 
an attractively-bound dimer in a combined periodic and 
parabolic potential with J/Q, = 140 and U = — 10 J. Initial 
state |^(0)} corresponds to the ground state Ix?) of the effec- 
tive Hamiltonian J5J with the dimer amplitudes af (0) shown 
in the left inset. Right inset is the projection £V I^l^)! 2 . 

and eigenstates for an effective dimer as 

Eg ra -2JW + 2\J JWM 2 ) (fc + i), (9a) 

| X f) « tf y £(2>>k\)- 1 '*e-&'H k (Z j )\lf), (9b) 

i 

where energies E® are relative to the dimer internal en- 
ergy (U - J< 2 >), = j$W/J(2) = jVWW k the 
discrete coordinate, and = cj |0) denotes a state 

with a single dimer at site j; obviously |1^) = |2j). The 
modified Bloch band -2 < Ej? < 2J^ for the dimer 
is restricted to the sites with 

W<fi. -^(l + = 1.3^, (10) 

thus containing N D = 2[jm ax J + 1 energy levels E® 
with < k < N D . The effective harmonic oscil- 
lator frequency at the bottom of the modified Bloch 
band is hcu D — 2V J( 2 )fK 2 ) and the dimer effective mass 
pP = h 2 /{2J^d 2 ) is large (j( 2 > < J) and positive. We 
have verified these conclusions by numerically solving the 
Schrodinger equation using the exact Hamiltonian (CQ) 
with the initial conditions corresponding to eigenstates 
(|9b[) of the effective Hamiltonian ([5]). As an example, in 
Fig. [4] we show the time evolution, or nearly complete 
absence thereof, of the system in the ground state of |8]), 

3 3 

(11) 

with energy Eg = -2 j( 2 ) + VT^H^T. 

We next turn to the case of strong repulsive interac- 
tion U > 0. The dimer tunneling constant is negative, 
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FIG. 5: Time evolution of density pj = (hj) for repulsively- 
bound dimcr in a combined periodic and parabolic potential 
with J/Q = 140 and U = 10 J. (a) Initial state |*(0)} is 
the ground state Ix?) °f attractive dimer, Eq. QT) , with 
amplitudes af(0) shown in the left inset, (b) Initial state 
|*l/(0)) is the ground state \xE) or repulsive dimer, Eq. (|14|l . 
with amplitudes af(0) shown in the left inset. Right insets 
are the projections Y^j |(2j|$)| 2 . 
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FIG. 6: Time evolution of (a) atom density pj = (hj), and 
(b) dimer density pf = (rhj ) ~ pj /2 in a combined periodic 
and parabolic potential with J/Q = 140 and U = 10J. Ini- 
tial state 1^(0)) corresponds to the dimer ground state \xx>) 
shifted by 3 sites from the trap center, (a) is the numerical 
solution of the Schrodinger equation with the exact Hamil- 
tonian {TJ, while (b) is obtained with the effective Hamilto- 
nian Inset in (a) is the projections ]TV |(2j |^}| 2 . 



J^ 2 ) < 0, corresponding to a negative effective mass fj, D 
[ill ]. As a result, \xq) in Eq. (fTTj) is no longer the ground 
state of Hamiltonian ([8]), as attested in Fig.[5ta). Rather, 
it is a highly excited state. To see this, consider for a 
moment a single particle in a flat lattice of N sites with 
J < 0. It follows from Eqs. (|2a|). (pb)) that the low- 
est energy state with Ef 1 _ 1 — — 2J cos[nN / (N + 1)] = 
-2|J|cos[7r/(7V + 1)] is 



N 

IXjv-i) = -AA^sin 



In 



N + l 



(12) 



Thus, in the limit of infinite lattice N — > oo, the 
ground state corresponds to the Bloch wave with quasi- 
momentum q = n. Returning back to the repulsively- 
bound dimer in the combined periodic and parabolic po- 
tential, we find that the low-energy eigenvalues are those 
of Eq. (|9aj) with the replacement — > \J^\, while the 
corresponding eigcnstates are given by 



lx?> 



i N Y j {2 k k\)- 1 ' 2 e-^V 2 H k {^)e^ \lf) 



(13) 



The ground state with Eg = -2|j( 2 )| + y/\j( 2 )\M 2 ) is 



then 



lx?> 



n( 2 ) 



7T 2 |J( 2 



E 1 



m\u\ 
^j 2 



^e-^-iy^j), (14) 



which is confirmed by our numerical simulations illus- 
trated in Fig.^b). Remarkably, the repulsive dimer ap- 
pears to be tighter bound than the attractive one. The 
symmetry between the cases of U < and U > is 
broken due to the presence of a parabolic potential. 

Finally, in Fig. [5] we show the dynamics of a dimcr 
wavepacket | , F), represented by the ground state Ix?) 
initially shifted by 3 sites from the trap center (for these 
parameters, jg ax ~ 4.9). Our simulations using the ex- 
act Hamiltonian (fTJ) and the effective Hamiltonian ([8]) 
yield practically identical results, which amount to peri- 
odic oscillations of the dimer wavepacket between the two 
sides of parabolic potential with period t d ~ 2tt/uj d = 
(Tth/2J)WU/D,. Numerical simulations for attractively 
bound dimers reveal similar behaviour but with con- 
siderably larger admixture of the single-particle states, 
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\ (lj\^)\ 2 < 0.2. This is another manifestation of the 
fact that the repulsive dimcr in a combined periodic and 
weak parabolic potential is bound tighter than the at- 
tractive dimer under the otherwise similar conditions. 



IV. CONCLUSIONS 

To summarize, in this paper we have studied coher- 
ent quantum dynamics of one and two bosonic parti- 
cles in a combined tight-binding periodic and shallow 
parabolic potential. Our studies are relevant to cur- 
rent experiments with cold alkali atoms in optical lat- 
tices and weak magnetic (or optical) traps 0]. After 
revisiting the single-particle problem, we considered ef- 
fective interaction-bound dimers recently realized in the 
experiment [l2j with strong repulsive atom-atom inter- 
actions. We examined both cases of repulsively-bound 
and attractively-bound dimers and identified similarities 



as well as marked differences in their static and dynamic 
properties. In particular, a rather counterintuitive fea- 
ture of the system revealed by the present work was that 
the repulsive dimers are bound stronger than the attrac- 
tive dimers, as far as their ground states and coherent dy- 
namics associated with low-energy states is concerned. In 
addition, we have shown that non-dispersive transport of 
carefully prepared atomic wavepackets can be achieved. 
As an extension of the present work, we plan to study 
dimer-monomer resonant collisions and entanglement of 
the resulting wavepackets. Our results may be pertinent 
to quantum information schemes with cold atoms in op- 
tical lattices. 
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